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In this article, we present a numerical study of the optical properties of a metal nanowire in¬ 
teracting with a localized cavity in a two dimensional photonic crystal. The nature of the modes, 
their wavelength and width are investigated as a function of the particle radius. It is shown in 
particular that a particle with diameter about the lattice constant presents very narrow resonances 
corresponding to hybrid photonic-plasmonic modes, where increased lifetime is attributed to the 
decrease of the radiative losses by interaction with the photonic crystal. These results open inter¬ 
esting applications in areas where narrow plasmonic resonances are required, as in LSP-resonance 
based biosensing. 


INTRODUCTION 


Localized surface plasmons (LSPs) are optical modes 
supported by metal nanoparticles, which result from the 
resonant coupling between the oscillating surface charge 
created along the particle’s surface by an incident electro¬ 
magnetic field, and the light wave that it scatters [T] . For 
a specific wavelength, depending on the particle shape, 
volume, material, as well as the optical characteristics 
of its environment, the extinction cross-section increases, 
and the intensity of light can be enhanced by several 
orders of magnitude in a length range around the parti¬ 
cle well below the wavelength of the incident light [2]. 
For all of these properties, LSPs are investigated due 
to both fundamental aspects and applications into nu¬ 
merous domains like biosensing EHS], surface enhanced 
Raman spectroscopy [6], photo-thermal therapy m , or 
metamaterials and optical cloaking [8] . In number of ap¬ 
plications, its is required to have a resonance as narrow as 
possible, for example in biosensing where the principle is 
to detect a small shift in the resonance wavelength when 
a target molecule is adsorbed on a metal nanoparticle 
surface: the resonance width determines the resolution 
of the sensor, and is often the main limitation in LSP- 
resonance-based sensing [9] . As the width of the plasmon 
resonance is determined by ohmic (which cannot be sup¬ 
pressed) and radiative losses, the narrowest resonances 
are obtained with the so-called ’’dark” plasmons modes: 
the distribution of electric charges on the particle surface 
is such as the particle does not possess any net dipole, 
and the radiative waves are cancelled out [T0l|TT] . Narrow 
plasmon resonances can be obtained as well by tailoring 
the environment of the particle, for example by making it 
interact with photonic crystals (PCs), as these structures 
present forbidden bands for radiative waves, which can 
be tuned by playing with their geometrical parameters. 
Recently, Solano and co-workers have shown [laiiH] that 
the absorbance of coated gold nanospheres and nanorods 
could be tailored when placed in a defect layer of a mul¬ 
tilayer planar resonator: enhancement of the absorption 
takes place within the band-gap at wavelengths close to 
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FIG. I: The investigated structure is a square photonic crystal 
with lattice parameter a, made of air-holes with radius r inside 
a Ti 02 matrix, where one of the hole is replaced by a metal 
nanowire of varying radius R. 


the one of the cavity defect, and narrowing of the res¬ 
onance results from the filtering of light by the cavity 
mode. Following that work, Wang et al m have shown 
theoretically that the same system could show greatly 
enhanced photo-induced change of the optical transmit¬ 
tance, based on the modification of the dielectric function 
of the gold nanoparticles on a subpicosecond time scale 
with an ultrashort light pulse. In that perspective, in¬ 
creasing the confinement of the localized plasmon mode 
of a nanoparticle using two-dimensional photonic crystal 
is of great interest for novel photonic devices. Hence, 
we propose in this article a detailed study of a 2D pho¬ 
tonic crystal interacting with an arbitrary size 2D metal 
nanoparticle. 

More precisely, we investigate in this article the opti¬ 
cal properties of a two-dimensional metal nanowire (NW) 
placed inside a localized cavity of a photonic crystal (Fig¬ 
ure [^. We aim at characterizing the optical modes of 
that system, by their nature, wavelength and resonance 
width, together with their electromagnetic field distribu¬ 
tion. The geometrical parameters of the photonic crystal 
have been chosen in order to have the localized plasmon 
wavelengths of the particle within the band-gap of the 
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PC. We expect that the radiative contribution to damp¬ 
ing rate will be strongly modified as the radiative waves 
emitted by the nanowire at the resonance are prevented 
from leaking inside the PC. The PC is made of a square 
grating of 2D circular air holes of radius r and lattice pa¬ 
rameter a, embedded inside a titanium dioxyde (Ti 02 ) 
matrix. This material has been chosen because of its 
transparency in the visible range of wavelengths, where 
its refractive index has been taken as constant and equal 
to 2.9 in the following. That high value allows ensuring 
a large index contrast between the material and the air 
holes, which allows obtaining large and deep band-gaps, 
that could not be reached with lower index materials like 
Si02. The metal 2D nanowire (NW) is a cylinder of ra¬ 
dius R made of gold, and is placed at the center of a 
localized cavity obtained by filling an air hole with Ti 02 . 
The moderate refractive index of Ti02 allows keeping the 
surface plasmon wavelengths of the nanowires in the vis¬ 
ible range. The data for the dielectric constant of gold 
are extracted from Johnson and Christy [15], and the cal¬ 
culations were performed using finite element methods 
(Comsol Multiphysics, Radio-Frequency module). In the 
first part, we present the properties of the gold nanowire 
embedded in an homogeneous Ti02 medium. In parallel, 
we show the band-gap tuning properties of the photonic 
crystal as a function of the grating filling factor, together 
with the calculation of the photonic modes sustained by 
the localized cavity. In the second part, we present how 
both the photonic modes of the PC and the particle’s 
plasmonic modes are affected when the metal particle is 
placed inside the cavity. The results are discussed as a 
function of the radius of the gold NW, first when it is 
much lower than the grating’s period, and second when 
both the radius and the period have similar values. 


OPTICAL PROPERTIES OF THE METAL 
NANOWIRE AND THE UNCOUPLED PC 

Au nanowire inside an homogeneous and infinite 
Ti02 medium 

Figure shows the extinction and absorption spectra 
of the gold nanowire in the homogeneous Ti 02 matrix 
as a function of its radius R. The particle have several 
high-order plasmon resonances which have their largest 
extinction signature between 600 and 700 nm. The reso¬ 
nances are labelled by their order m, m = 1 correspond¬ 
ing to the dipolar mode, m = 2 to the quadrupolar mode 
and so on. As the refractive index of the embedding 
medium is quite high compared to more usual materials 
in plasmonics, like silica or water, the LSP modes sup¬ 
ported by the particle can have a large order m, up to 
m = 7 for the 150 nm radius particle. The quadrupo¬ 
lar mode appears in the visible range of wavelengths for 
a diameter of about 40 nm, for which the dipolar reso¬ 


nance is red shifted and barely noticeable. The absorp¬ 
tion map (Fig.j^b)) is quite different from the extinction 
map (Fig. [^a)), due to from the fact that the contribu¬ 
tion of the scattering to the extinction increases with the 
particle size. As can be seen on Fig. the absorption 
is comparable to the scattering for a small particle of ra¬ 
dius about 10 nm, but, as the particle is rapidly large 
compared to the wavelength of the light in the medium 
(far from the electrostatic approximation), the scattering 
largely dominates when the particle diameter increases, 
as it can be verified for R = 85 nm. This is the reason 
why in Fig. I^the extinction is high for larger radii, while 
the absorption appears low. 

(a) Extinction 
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FIG. 2: Extinction (a) and absorption (b) plotted as a func¬ 
tion of the wavelength and the nanowire’s radius (arbitrary 
units, maximum in red, minimum in blue). The labels m 
indicate the order of the different localized plasmon mode. 


Photonic crystal 

The geometrical parameters of the photonic crystal 
have been chosen in order to tune its absolute band-gap 
to the 600 nm / 700 nm wavelength domain, where the 
nanowire’s LSP modes are the most efficiently excited. 
As the refractive index is fixed, the only two degrees of 
freedom are the period a of the grating and the radius 
r of the air hole, which both allow to change the cen¬ 
tral wavelength, the width and the depth of the band 
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FIG. 3: Comparison between the extinction (black), absorp¬ 
tion (red) and scattering (blue) cross-sections for a particle 
of radius = 10 nm (dashed line) and a particle of radius 
R — Sb nm (solid line). 


gap. We have fixed the period to a = 320 nm, while 
figure [^a) shows the evolution with r/a of the transmis¬ 
sion spectra of a planewave in normal incidence onto one 
of the (10) direction of the grating, limited to 7 lines of 
holes. The polarization of the incident wave is TM, the 
incident magnetic field being along the invariance direc¬ 
tion. In this transmission calculations and the follow¬ 
ing, periodic boundary conditions have been applied in 
the x-direction, while perfectly matched layers have been 
applied in the incidence and emergence media, perpen¬ 
dicular to the incident planewave direction. The trans¬ 
mission spectrum has been computed by integrating the 
Poynting vector flow through a line parallel to the PC 
in the emergence medium. It appears that the structure 
supports several band gaps, and that for R/a ^ 0.42 a 
large band gap overlaps with the target wavelength do¬ 
main (blue area). In order to check the consistency of 
the transmission calculations, we have compared these 
results with a full modal dispersion calculation, as shown 
on Fig. [^b), for a = 320 nm and r = 0.42a. We can see 
that a complete band gap exists in the reduced frequency 
domain a/X G [0.4524; 0.4841] (A G [661 nm; 707 nm] 
with a = 320 nm), whereas a partial band gap in the 
TX direction exists in the interval a/X G [0.452; 0.554] 
(A G [578 nm; 707 nm] with a = 320 nm). 

A cavity is then created by filling one of the air holes 
with Ti02. On Figure [^a), the dispersion calculation on 
a 7x7 supercell shows that two localized photonic modes, 
respectively (a) and (s), are found inside the complete 
band-gap for reduced frequencies a/X respectively equal 
to 0.48 (A = 668 nm) and 0.46 (A = 694 nm). The 
distribution of the electric field amplitude is shown on 
the colormaps. Fig. [^a), right. The lowest wavelength 
mode (a) has an annular shape, its field lines form closed 
loops around the center of the cavity where the inten¬ 
sity reaches a minimum. We can notice that the mode 
is antisymmetric with respect to the two principal axes 
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FIG. 4: (a) Transmission of a TM-planewave in normal inci¬ 
dence on a PG made of seven lines of air holes with radius r 
and period a = 320 nm, as a function of r/a. Maximum trans¬ 
mission is white, minimum is black; (b) Band diagram of the 
PG along the high symmetry axes FA, AM and MF of the 
first Brillouin zone, with period a = 320 nm and r/a — 0.42. 
The gray areas represent the absolute band gap whereas the 
red areas represent the band gap in the FA direction. 


of the photonic crystal cavity. The highest wavelength 
mode (s) has a very different distribution of electric field, 
and shows several local maxima inside the cavity. It is 
symmetric with respect to the main axes of the photonic 
crystal cavity and reaches as well a minimum amplitude 
at the center of the cavity. Figure [^b) shows the trans¬ 
mission spectra of a TM-planewave normally incident on 
the (10) direction of the grating, still limited to 7 rows. 
The wavelength and electric field distribution of the main 
resonance is consistent with the ring mode (a) found with 
the dispersion calculation. The mode at 694 nm is not 
excited, because, as the incident light is normally inci¬ 
dent onto the photonic crystal, it is antisymmetric with 
respect to the vertical direction, and cannot couple to 
the 694-nm-symmetric mode. However, the same cal¬ 
culation performed with a small incidence angle breaks 
the symmetry of the incident field and allows to recover 
the corresponding transmission peak (not shown) at the 
wavelength predicted by the dispersion calculation. 
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FIG. 5: (a) Band diagram of the 7x7 supercell of the same 
PC with a localized cavity. The wavelength and electric held 
distributions have been indicated for a = 320 nm; (b) Trans¬ 
mission of a TM-planewave in normal incidence on the same 
PC, limited to 7 lines along the direction of the incident wave. 


COUPLED SYSTEMS 


In this section, we investigate the optical behaviour of 
the gold nanowire placed at the center of the localized 
PC cavity. The absorption spectra have been calculated 
as a function of the particle radius, by integration of the 
power losses on the cross section of the nanowire, and 
normalization to the incident energy flow on the 7x7 su¬ 
percell. Hence, reflection (i7), transmission (T) and ab¬ 
sorption {A) obeys the energy conservation i7+T-j-H = 1. 
Two very different behaviours have been obtained, ac¬ 
cording to the radius of the particle compared to the 
grating’s period. 



A (nm) A (nm) 


FIG. 6: Transmission and absorption spectra of the system 
composed of the gold NW of radius R placed at the cen¬ 
ter of the cavity, for wavelengths close to the cavity ring¬ 
mode position. For readibility purpose, spectra have been 
shifted by 0.005 per radius increment. Maxima of transmis- 
sion(absorption) are: 8.5(0.24)% (R = 10 nm), 6.0(2.0) % (R 
= 20 nm), 2.5(2.6) % (R = 30 nm), and 0.9(1.5) % (R = 40 
nm). 


Case where R « a 

We consider first the case where the particle has a ra¬ 
dius R small compared to the lattice parameter a. Fig¬ 
ure shows the evolution of the transmission peak of 
the ring mode (a), from radius 10 to 50 nm. The wave¬ 
length of the mode is blue-shifted when R increases, while 
its transmission rapidly decreases. For comparison, the 
absorption spectra have been added, which shows that 
both absorption and transmission reach comparable low 
values, and that most of the energy is then reflected. Ta¬ 
ble [T] shows the different resonant wavelengths, together 
with the electric and magnetic field intensity distribution 
for R < 40 nm. It appears that the particular ring struc¬ 
ture of the unperturbed PC cavity modes persists for val¬ 
ues of R up to few nanometers. Indeed, as the intensity 
of the PC-cavity ring mode is almost zero at center, its 
shape and wavelength is unchanged in the limit of very 
small particle {R = 10 nm). However, when R increases, 
we can intuitively understand the decrease of the mode 
wavelength by the fact that the light is pushed away from 
the area occupied by the metal particle: this leads to a 
decrease of the area accessible to the photonic mode, and 
then a shortening of the resonance wavelength. 


Moreover, if the shape of the mode is unchanged for small 
radius, it is not the case for a particle radius larger than 
50 nm, where the light ring is much closer to the adja¬ 
cent air holes. This breaks the revolution symmetry of 
the ring mode, and for R > bO nm the ring mode does 
not exists anymore. Finally, let us emphasize that this 
mode, when the particle radius is much smaller than the 
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( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 


R(nni) 

0 

10 

20 

30 

40 



A(nm) 

668.6 

668.5 

667.6 
664.5 
658.9 


TABLE I: Summary of the wavelengths and field distributions 
for the different peaks labelled (1) to (5) on Fig.[^ 


lattice constant, is purely photonics, without plasmonic 
characteristics, as the electric field is, in every point of 
the metal interface, parallel to it, and then cannot create 
the surface charge necessary to excite surface plasmon 
waves. 


Case where 2R is comparable to a 

The opposite regime is when the particle diameter is 
comparable to the lattice parameter. In that case, the 
light is squeezed in the very narrow space in between the 
gold interface and the air holes, and strong coupling oc¬ 
curs between the gold nanowire and the photonic crystal. 
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FIG. 7: Transmission and absorption spectra of the coupled 
system with the NW’s radius R comparable to the grat¬ 
ing’s period a. For readability purpose, spectra have been 
shifted by 0.001 per radius increment for transmission and 
0.5 per radius increment for absorption. Maxima of transmis- 
sion(absorption) are: 0.03(3.0)% (R = 100 nm), 0.33(9.5)% 
(R = no nm), 0.75(13.9) % (R = 120 nm), 0.60(13.6) % (R 
= 130 nm), 0.34(9.7) % (R = 140 nm) and 0.15(6.1) % (R = 
150 nm). 

Figure shows the evolution of the transmission and 
absorption spectra for particle radius varying from 100 
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TABLE 11: Structure of the compounds modes (1) and (2) for 
R = 150 nm. 


nm to 150 nm. Maximum values, indicated in the cap¬ 
tion, show that the incident energy is mostly reflected, 
but a significant part is absorbed (up to 14% for R = 
120 nm) whereas less than one percent is transmitted. 
This is a very different behaviour than for small radius, 
where both maximal absorption and transmission are of 
the same order of magnitude. For R between 100 and 
130 nm, only one mode is excited between 600 and 700 
nm, labelled (1). It is red-shifted with R until R = 110 
nm, while for larger radius it blue shifts. Similarly, a 
second peak, labelled (2), appears for R = 140 nm at 
a wavelength of about 660 nm, and is blue-shifted with 
R. This blue-shift of these two modes with large values 
of R might be surprising as the localized plasmon modes 
of a single nanowire usually shift to the red with diame¬ 
ter (see Fig. [^. This is actually what is observed when 
R <110 nm. when R becomes larger, the reversed evolu¬ 
tion of the wavelength can be explained by the decrease 
of the effective refractive index surrounding the particle: 
for wider NW, the air holes are closer to the wire’s edges 
and the portion of low refractive index (air holes) in the 
area occupied by the optical mode is larger. In that case, 
the blue shift due to the lowering refractive index is larger 
than the red shift due to the increasing diameter. 


The distribution of the electric and magnetic field am¬ 
plitude, together with surface charges (plotted at the 
time of maximal amplitude) are displayed on Table 
for the modes (1) and (2) and R = 150 nm. In both 
cases, the distribution of the electric field looks complex 
and difficult to interpret. However, the magnetic field 
structure is much simpler and shows clearly a similarity 
between the mode (1) and the octupolar mode (m = 4) 
of the single nanowire in the homogeneous Ti 02 matrix: 
the main qualitative difference is at the proximity of the 
air holes where one intensity maxima out of two around 
the particle is lower than its nearest neighbors. The simi¬ 
larity with the octupolar mode is confirmed by the repar¬ 
tition of the surface charges which clearly demonstrates 
that this mode is antisymmetric with respect to the two 
principal axis of the photonic crystal. 
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The mode (2) has a very different structure as indi¬ 
cated by the magnetic field distribution, as the maxima 
around the nanowire surface are evenly distributed, ex¬ 
cept on the axis perpendicular to the excitation direction 
where the light intensity is zero. This can be under¬ 
stood by the surface charges distribution which is anti¬ 
symmetric compared to the vertical main axis whereas 
it is symmetric compared to the horizontal axis. This 
results in a zero of field along the nanowire’s horizontal 
diameter. This last mode is clearly different from the sin¬ 
gle nanowire modes and as a direct consequence results 
intrinsically from the coupling of the localized plasmon 
of the particle with the photonic crystal. As both modes 
(1) and (2) result from the strong coupling of the NW 
and the PC, we will call them PC-plasmon modes in the 
following. 



600 620 640 660 680 700 

A (nm) 

FIG. 8: Comparison between the absorption spectra for the 
single nanowire in infinite Ti 02 medium (a) and the nanowire 
of radius 140 nm in the PC cavity. 

We compare on Fig. |^the modes (1) and (2) for the 
R = 140-nm particle (a) with the modes of the isolated 
nanowire with the same radius in the homogeneous 
matrix (b). We can clearly see that the PC-plasmon 
modes are much narrower. More precisely, mode (1) and 
(2) have a similar full width at half maximum (FWHM) 
of about 6 nm, which corresponds to a quality factor of 
113 for mode (1), whereas the m = 6 mode at 680 nm 
for the isolated nanowire has a FWHM of 36 nm with a 
quality factor of 19. This reduction of the mode width 
by a factor of about 6 is attributed to the lowering of 
the radiative contribution of the electromagnetic field 
scattered by the nanowire at the resonance, which is 
forbidden to propagate into the photonic crystal. This 
phenomenon should be responsible for the decreases of 
the radiative losses, making the plasmon mode width 
limited principally by its ohmic losses. 


Finally, the last figure summarizes the different phe¬ 
nomena involved in the coupling of the metal particle 
with the photonic crystal for a radius varying between 0 
and 150 nm by step of 10 nm, with few representative 
magnetic field distributions. Transmission and absorp¬ 
tion spectra have been plotted side-by-side in log scale. 
As previously explained, for small radius, the ring mode 
is not strongly affected by the particle and shifts to the 
blue {R branch), whereas the PC-plasmon modes (1) and 
(2) appears for radius larger than 110 nm. In between, 
we can see one mode which is linearly red-shifted with 
R: it corresponds to the octupolar, m = 4, mode of 
the nanowire in the homogeneous matrix. The hexapo- 
lar m = 3 mode appears at A = 680 nm for R = 60 nm, 
but shifts beyond 700 nm for larger radius. Actually, 
for small radius, these modes are not strongly affected 
by the photonic crystal as the plasmon mode does not 
strongly overlap with the air holes of the photonic crys¬ 
tal. A transition occurs at about R = 100 nm where these 
modes start to combine with the photonic crystal to give 
the PC-plasmon modes (1) and (2) previously described, 
where they show much narrower than modes m = 3 and 
m = 4. 



A (nm) A(nm) 

FIG. 9: Evolution of absorption in the metal nanowire, as a 
function of the wavelength and the particle radius. 


SUMMARY 

In this work, we have investigated the optical prop¬ 
erties of a 2D metal nanowire placed at the center of a 
localized cavity in a 2D photonic crystal. The nature of 
the modes, their wavelength, together with their evolu¬ 
tion with the particle radius is very dependent on how 
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the diameter compares to the grating period. The ring 
photonic mode of the PC cavity has his shape essentially 
unaffected by a particle of few ten nanometers diameter, 
while its wavelength blue-shifts. However, for a parti¬ 
cle with diameter about the lattice constant, the coupled 
system present very narrow resonances corresponding to 
compound PC-plasmon modes, where the narrowing of 
the resonance is attributed to the decay of the radia¬ 
tive losses of the plasmon particle by the photonic crys¬ 
tal. These results open interesting applications in ar¬ 
eas where tailoring the plasmon response of nanoparti¬ 
cle is required, as in LSP resonance biosensing or high¬ 
speed photonic devices. Beside, such a system is a choice 
platform for investigating the coupling of localized plas¬ 
mon modes with elastic waves, as it has previously been 
demonstrated that PC can leads to enhanced interaction 
between phononic and photonic waves [min]. Similarly, 
such a compound photonic crystal/metal particle sys¬ 
tem could show enhanced interaction between elastic and 
plasmonic waves p!8f[2Q] . allowing to achieve mechani¬ 
cal control of the plasmon response at a sub-micrometer 
scale. 
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